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FUNCTIONS OF LINES IN BALLISTICS* 



GILBERT AMES BLISS 

In the theory of ballistics it is customary first of all to set up the differ- 
ential equations describing the motion of a projectile under the assumptions 
that no wind is blowing and that the densities of the air at different altitudes 
are normal. If the initial velocity and the initial direction of flight of the 
projectile are given, it is then possible to compute the coordinates of the 
projectile as functions of the time t satisfying .these differential equations. 
The path of the projectile is different, however, when a wind is blowing or 
when the air densities are not normal, and the form of the differential equations 
themselves must be changed to account for such disturbed conditions. The 
new differential equations involve the velocity w ( y ) of the wind in the direc- 
tion of the plane of fire at the altitude y , and the correction k ( y ) for the 
density at the altitude y . The range X thus becomes a function X [ w ( y ) , 
k ( y ) ] of the functions w(y) , n(y) , or in other words a function of lines, 
and the first differential of X in the sense of the theory of functions of lines 
is the first order effect of wind and density upon the range. 

The range is the quantity of greatest interest to the ballistician in most 
problems of artillery fire. But for anti-aircraft guns one must also investigate 
the corrections to the three coordinates x , y , z of the projectile at the ends of 
successive short intervals of time t along the trajectory. It is clear that 
these coordinates will also be functions of lines in the sense described above 
for the range X . 

In the following pages a method for computing the differentials of such 
functions of lines is described. The formulas found are the results of appli- 
cations, to differential equations of the forms used in ballistics, of more general 
results deduced by the writer in a preceding paper, f The proofs here given, 
however, have been made independent as far as possible of those for the 
general case. The methods used have proved to be of value in the practical 
computations of the differential corrections hitherto deemed necessary for 
trajectories, and as the applications of the theory of ballistics become more 

* Presented to the Society, September, 1919. 

t Differential equations containing arbitrary Junctions, these Transactions, vol. 21 
(1920), p. 79. 
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refined the formulas will very likely be of value in the computation of other 
corrections also. The writer has already published accounts of the method 
with such mathematical details as are necessary for the practical ballistician.* 
The purpose of the present paper is to show the connection between the 
ballistic theory and the theory of functions of lines. The method given for 
finding the first differential of a function of lines is one which could be applied 
to many cases of such functions defined by differential equations. 

1. The fundamental formula for differential corrections 
Consider a set of differential equations of the form 
x" =f(t, x, y, z, x' , y' , z') , 

(1) y" = g(t, x, y, z, x' , y' ,z'), 
z" = h(t,x,y,z,x',y' ,z'), 

and having an initial solution 

(C) x = x(t), y = y(t), z = z(t) (* x <*<*,) 

for which the functions x(t) , y (t) , z(t) have continuous derivatives of the 
first and second orders on the interval t\ t 2 , and which passes through an 
initial element ( t , x , y<> , %o , x'o , y'a , z'o ) • In a neighborhood N of the values 
(t, x, y , z, x' , y' , z') on this solution the functions/, g, h are supposed to 
have continuous partial derivatives up to and including those of the second 
order at least. 

If the functions/, g , h are replaced by three others /, g , h, having similar 
properties in the neighborhood N , the solutions of the equations 

x" = f(t, x, y, z, x' , y' , z'), 

(2) y" = g(t, x, y, z, x' , y' , z'), 
z" = h(t,x, y, z, x' , y' , z'), 

will be uniquely determined by an initial element (to, x , yo, %o, x a , y , z'o) , 
and if U,xa,yo, Zo are kept fixed x may be regarded as a function x [ x , y'o , z'o, 
/, g , h] with similar expressions for v and z . It is natural to expect that by 
taking the elements of the set [ x , y , z'o, f , g , h] sufficiently near to those 
of the set [ x' , y' , z'<, , f , g , h] belonging to the solution C the three functions 

* A method of computing differential corrections for a trajectory, Journal of the 
United States Artillery, vol. 51 (1919), p. 445; The use of adjoint systems in 
the problem of differential corrections for trajectories, ibid., p. 296; Differential corrections for 
anti-aircraft guns, (not yet published). 



1920] FUNCTIONS OF LINES IN BALLISTICS 95 

x , y , z may be made to differ as little as is desired from the functions x ( t ) , 
y(t),z(t) defining that solution. That this is the case follows readily from a 
theorem of the first paper by the writer referred to in the introduction above,* 
provided that the equations ( ) are written in the normal form 

dx _ , dy _ / dz f dx' dy' dz' 

dt =x ' dl = y ' ~dt = z ' ~dT = *' ~dT = g ' ~dl = h - 

If Aa; denotes the difference x — x , with similar notations for the other vari- 
ables, the quotients Ax/mi , Ax' '/mi , and the similar ones for y and z , are all 
bounded as mi approaches zero, m x being the maximum of | Axo | , |Ay u |, 
|A»;|, |A/|, \Ag\, |AA|.f 
By subtracting equations (1) from equations (2) a relation of the form 

Ax" = f(t,x,y,z,x',y',z')-f(t,x,y,z,x',y',z') 

+ Af(t, x, y, z, x' , y' , z') 

is found for Ax, and similar ones hold for Ay and Az. This equation may 
also be written in the form 

(3) Ax" = f x Ax + f v Ay + /, Az + /,' Ax' + f y > W + /.' Az' + Af + <j> , 

where 

<t> = f(t, x, y, z, x' , y' , z') - / - f x Ax - f v Ay -f z Az 

- W A*' - W A V' ~ U' As' + Af(t,x,y,z,x',y',z')- Af. 

Similar equations hold for Ay and Az with/, <j> replaced by g, \p and h, x- 
It is understood that arguments not indicated in /, g, h and their derivatives 
are always those belonging to the curve C. By applications of Taylor's 
formula with remainder terms of the second degreej to the first eight terms 
of the expression for <f>, and also to the last two, it is clear that <j> is of the 
second degree in Ax, Ay , Az, Ax' , Ay' , Az' and the derivatives Af x , Af y , Af z , 
Af x ' , Af u ' , Afz' . If m 2 is the maximum of mi and | Af x \ , \ Af y \ , \ Af 2 \ , \ Af x ' | , 
| Afy' | , | Afz' | , then <f>/rnl is bounded as rtii approaches zero, and <j> = m| ( <j>/m\ ) 
itself may therefore be regarded as being of the second degree in m 2 , as will 
also be the case for \p and x • 

Equation (3) and the similar ones for Ay" and Az" may be written in the 
normal form 



* Loc. cit., p. 85. 
t Loc. cit., p. 89. 

+ The formula with integral form of the remainder terms is the most satisfactory. See 
Jordan, Cours d'analyse, vol. 1, 2d edition, p. 247. 
Trans. Am. Math. Soc. 7 
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dt 



Ax = 



Ax' , 



d A 
dt Ay = 

d A 
-dl AZ = 



At/', 



Az', 



(4) 



d_ 

dt 

<L 

dt 

d_ 

dt 



Ax' = f x Ax + f y Ay + U Az + f x > Ax' + f v > Ay' + f t > Az' + A/ + <j> , 
Ay' = g x Ax + g y Ay + g* Az + g x > Ax' + g y > Ay' + g z > Az' + Ag + \j/ , 
Az' = h x Ax + h v Ay + h z Az + h x ' Ax' + h v ' Ay' + h z ' Az' + Ah + x ■ 



The system adjoint to this is by definition the system 



(5) 



-X' = 


fz P + g x cr + h x r , 


— p' = 


fy P + 9v a + Kt, 


— v' = 


fz P + g z a + Kt , 


-p'=X 


+ fx'p + g x ' o" + h x ' t , 


— <x' = p 


+ fy'P + 9v'a + K' T > 


— t' = 


V +/«' P + 9z'<7 + h z ' T, 



formed by using the columns of coefficients of the system (4), with signs 
changed, in place of the rows. If X, p., v , p , a , t is a solution of the last 
system, then by multiplying the equations of the system (4) by them, re- 
spectively, and using equations (5), it follows that 



dt 



(\Ax + nAy + vAz + pAx' + a Ay' + tAz' ) 



= pAf + aAg + rAh + p<j> + a\p + t\ 
Furthermore by integrating this last result from t to t , the relation 
XAa; + pAy + vAz + pAx' + a Ay' + tAz' 

= [XAa; + pAy + vAz + pAx' + a Ay' + tAz'] (=(o 



(6) 

is found 






+ (pAf + aAg + rAh)dt + I (p</> + <rf + T\)dt 



f 
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In his last equation the second integral is of the second degree in the 
maximum ra 2 , since <f>, \f/ , x have this property, as explained above. The 
first two terms on the right may therefore be regarded as the differential of 
the sum on the left with respect to the differences Ax' , Ay' , Az' , Af,Ag,Ah. 
If t is considered as fixed for the moment, and if (X, p., v, p, a, t) are so 
chosen that their initial values at the particular value t are (1,0,0,0,0,0), 
then the formula (6) gives an expression for Ax in which the first order terms 
are the first differential 6a; of the function x in terms of the differences Axi , 
A^o, Az'o, Af , Ag , Ah. In a similar manner the differentials of y and z are 
determined. It is justifiable then to deduce a differential formula from (6) 
which takes the form 

X5x + tidy + vbz + phx' + <r§y' + tSz' 

(7) # r 

= [ pAx'o + <rAy'» + tAz ]'='» + I ( pAf + aAg + rAh ) dt . 

It is the equation which expresses the equality of the first order terms in 
equation (6) . The existence of the differentials of x , y , z , and their derivatives 
is also an immediate consequence of the theorems of the paper referred to 
above.* 

The last equation is fundamental for what follows. It gives the first order 
terms with respect to the variables Ax , Ay' , Az'o, Af, Ag , Ah of every ex- 
pression of the form 

(8) XAz + p.Ay + vAz + pAx' + a Ay' + rAz' 

for which (X,ju, v , p , a , t) is a solution of the adjoint system (5), and it 
may be called the first differential of such an expression with respect to the 
variables x' a , y' , z' and the functions f,g,h. If the first order terms are 
desired at a fixed value of t for an expression of the type (8) with arbitrarily 
chosen coefficients, it is only necessary to insert in the equation (7) a set of 
solutions \,p,v,p,(r,T of equations (5) with initial values at t equal to the 
coefficients in the expression (8). On the other hand, if t is variable, then 
six linearly independent solutions of equations (5) give six independent equa- 
tions of the type (7) from which the six differentials dx , dy , 5z, 5x' , by' , 8z' 
can be calculated. It is to be explained in the following sections how these 
considerations can be applied to calculate the differential corrections for 
trajectories. 

2. Formulas for range corrections and deflections in ballistics 

The differential equations of the trajectory of a projectile have the 
form 



* Loc. cit., p. 89. 
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X" = - R- = - Fx', 

V 

(9) V"= -R^-g= -Fy'-g, 

z" = - R- = - Fz', 
v 

where R is the retardation of the projectile in 
the direction of the tangent due to the resis- 
tance of the air, g is the gravitational con- 
stant, and F merely a convenient notation for 
R/v . The value of R used by Moulton is 

R = vF{v,y) = v — ^ , v = W + y + z , 

where vG(v) is the resistance of the air to a standard projectile moving at 
sea level with velocity v , C is the ballistic coefficient which enables the formula 
to be fitted to projectiles of other types, and e~ av is the altitude function 
which accounts for the diminution in the retardation as the altitude y increases 
and the density of the air diminishes. The values of G ( v ) have been tabu- 
lated for different velocities from observations of actual firings and the table 
has been smoothed by means of an analytic function which agrees with the 
observed data with sufficient accuracy. The value of the constant a = .000106 
for use with the meter-second system of units, is the result of meteorological 
observations extending over many years. It is so chosen that the function 
e~ av represents the ratio of the density of the air at the altitude y to the density 
of the air at sea level under so-called normal or average conditions. These 
normal conditions rarely exist at an actual time of firing, and the variations 
from them must be taken into account as will be explained in the following 
paragraph. 

It is presupposed that a trajectory 

(10) x = x(t), y = y(t), z = (O^t^T) 

has been computed with the help of equations (9), having given initial values 
x = y = z = , x' = x' , y' = y'a, z' = Oo at the time t = . The values 
of z will all be zero since the equations (9) have one and but one solution for 
the prescribed initial values, and since such a solution can be found by setting 
z = in the equations and integrating the first two for x and y . The time of 
flight of the projectile will be denoted by T and the range by X . The actual 
computation may be most economically effected in a particular case by the 
methods of approximation devised by F. R. Moulton and his associates in 
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the ballistic work at Washington,* the result being a table of values of the 
functions (10) at the ends of short intervals of time from t = to t = T . 
The problem of the present paper is then the determination of the corrections 
which must be applied to the values of these functions, and in particular to 
the range X , in order to account for disturbances of the trajectory of various 
sorts, some of which will be explained below. 

If the second members of equations (9) are denoted by/, g , h, and if they 
are replaced in any way by new functions 

/=/ + A/, g = g + Ag, h = h + Ah, 

in order to account for disturbances of the trajectory, then the results described 
in Section 1 will be applicable. The linear combinations (7) of greatest 
interest to the ballistician are the deflection Sz ( T ) and the differential SX 
of the range X . To find a first expression for the latter one may proceed 
intuitively as follows. The coordinates of the projectile on the disturbed 
trajectory at the time T are ( X + Ax , Ay , Az ) and the value 

AX = [Ax + cotuAy]'= T , 

where co is the positive angle of fall, 
is approximately correct for AX, as 
one sees readily from the figure which 
shows the projections of the original 
and the disturbed trajectory on the 
plane of fire. The value of the dif- 
ferential SX is found by replacing the 
differences Ax and Ay by differentials 
in this result. 

A more satisfactory method of determining SX from the standpoint of the 
mathematical analyst is, however, to apply the theorems on implicit functions 
of the paper referred to above.f The equation 

y[ T, x'o, y'o, z ,f, g, h] = 0, 

where y is the function of Section 1 above for the special case of equations (9), 
has an initial solution [ T , x' a , y' 0> 0,/, g, h] at the end of the trajectory (10), 
in a neighborhood of which the implicit function theorem is applicable. The 
solution T [ x'o , y' , Zo , f , g , h] of this equation is the time of flight on the 
disturbed trajectory and the expression 

X[x , y'o, lo,/, g, h] = x[ T , x'o, y'o, z ,/, g, h] 

* For a sketch of his methods see Journal of the United States Artillery, 
vol. 51 (1919), p. 40. 
t Loc. cit., p. 90. 
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is the corresponding range. According to the results of the paper just cited 
the differential of X is 

(11) SX = \ 5x ~-,8y 1 = [8x + cot o8y] t=T , 

where Sz and 8y are the differentials of x and y with respect to x , y' , z , 

f, g, h. 

The differential equations corresponding to the system (5) for the special 
case of this section are readily found to be 



(12) 



= 0, 

= — ax' Fp — ay' Fa , 

= 0, 



= -n + x'y'F^p+F^l+y'^y, 
= - v + Ft. 



To find 8X from formulas (7) and (11) a solution {\, p., v , p , a , t) oi this 
system with initial values ( 1 , cot a> , , , , ) at t = T is needed. For 
such a solution it is clear that v = r = so that the equation (7) gives 



(13) 



8X = p(0)Ax' (> + <T(0)Ayo+ ( pAf + <rAg ) dt . 

Jo 



Similarly to find 5z ( T) a solution of equations (12) must be used with initial 
values (0,0, 1,0,0,0) at t = T. The functions X , ju , p , a are all identically 
zero in this case since the first, second, fourth, and fifth of equations (12) are 
linear and homogeneous in X , p. , p , a alone and have only the solutions X = p. 
= p = a = with the initial values ( , , , ) . Formula (7) then gives 



(14) 



5z ( 



T) = r(0)Az' + rAhdt. 
Jo 



The function r in the last formula can be expressed in terms of the known 
function x(t) of the trajectorv (10). For F = - x"/x', from the first of 
the equations (9), and consequently the last of equations (12) gives 

t'x' + rx" = - x' , 

since now v = 1 . Integration of this result, and evaluation of the constant 
of integration by setting t = T , gives 

X -x(t) 



(15) 



x'(t) 
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The expressions (13) and (14) are the first differentials of the range X and the 
deflection z(T) thought of as functions of the ordinary variables x' , y' , z' and 
of the functions f , g, h. They give respectively the first order range correction, 
and the first order deflection, corresponding to a disturbance of the trajectory 
which can be accounted for by altering the initial conditions, or by changing the 
functions f , g , h in the second members of equations (9) into f + Af, g + Ag, 
h + Ah, or both. 

3. Applications to special corrections 

The formulas (13) and (14) are applicable to a very great variety of prob- 
lems. Suppose, for example, that experiments had been devised which 
gave a new and more accurate table of values of the function G ( v ) than the 
one at present used. It would not be necessary to re-compute completely 
the ranges corresponding under normal conditions to given elevations. For 
the new value of G could be inserted in equations (9) and the corrections 
A/, Ag, Ah computed. The formula (13) would then give first order cor- 
rections which could be applied to revise the ranges so that they would be in 
accord with the new retardation law. 

The most important applications of the formulas are, however, in computing 
the differential corrections due to variations in the initial conditions, variations 
from normal in the density of the air, or wind. The direction and velocity 
of the wind at different altitudes is determined in practice by observing the 
motion of a small balloon through two surveying instruments at the ends of a 
measured base line, and its components at the altitude y in the plane of fire 
and at right angles to it will be denoted hereafter by w x ( y ) and w t (y) . The 
variations from normal in the density of the air at different altitudes is deter- 
mined by means of airplane observations and may be accounted for in equa- 
tions (9) by using an altitude function of the form e~ ay [ 1 + k ( y ) ] instead 
of e~ ay . 

When a wind is blowing the retardation of the projectile is no longer in 
the line of the tangent to the trajectory, but takes place in the direction of 
the motion of the projectile with respect to the air, and is determined in 
magnitude by the velocity V of the projectile with respect to the air. This 
velocity is the vector sum of the velocity of the earth with respect to the pro- 
jectile and of the velocity of the air with respect to the earth, and its com- 
ponents and magnitude are 

- x' + w x (y), - y' , - z' + w z (y), 



U = V[x' - w x (y)f + y' 2 + [*' - w z {y)f, 

where x, y, z represent the coordinates of the disturbed projectile. The 
differential equations which determine the trajectory disturbed by wind or 
abnormal density are now 
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x" = -F[x'-w x (y)], 

(16) y"=-Fy'-g, 

z" = -F[z' -w,(y)], 
where 

(17) F = ±G(V)e-"i[l + K(y)]. 

The solutions of these equations are functions x , y , z of the arguments t , x' , 
y' u , z'u, w x (y) , w z (y) , ic(y) , and the first order deflection 8z ( T ) and range 
correction 8X are given by formulas (14) and (13) when Af ,Ag, Ah are replaced 
by the first order terms of their expansions in powers of w x ( y ) , w z ( y ) , n{y) 
which are described below. 

The second members of equations (16) are the functions previously desig- 
nated by /, g , h, and the differences A/ = / — /, Ag , Ah can be readily calcu- 
lated. In the formula (7) their arguments are the elements t,x,y ,z,x' ,y' ,z' 
of the known solution C , and in the formulas (13) and (14) it is to be re- 
membered that they must similarly be the analogous elements of the solution 
(10) already computed. It follows then that to first order terms in w x (y), 
w*{y), n(y) 

Af = F^l + x' i ~ S jw x (y) -x'F K (y) + ■■■, 

C 
Ag = x' y' F ^w x (y) - y' F n(y) + 

Ah = Fw z {y) + ■■■. 

When these values are substituted in (13) and (14), and use is made of the 
equations (12) and (15), it is found that 

(18) SX = p(0)Ax' o +<T(0)Ay' + f (p' + l)w x (y)dt + - f »' K (y)dt, 

Jo a Jo 

(19) Sz(T) =-tAzo+ f (r'+l)w z (y)dt. 

%o Jo 

The expressions (18) and (19) are the first differentials of the range X and 
the deflection z(T) thought of as functions of the ordinary variables x' a , y' 0> z'o 
and of the wind- and density-functions w x (y) , w z (y) , n(y) . They give, 
respectively, the first order range correction and the deflection due (I) to a variation 
in the initial velocity with projections Ax'o, Ay' a , Az'o on the coordinate axes, 
(2) to a wind at the altitude y with components w x (y) , w z (y) in the plane of 
fire and at right angles to it, and (3) to a variation from normal in the density 
ratio of the air which at altitude y is equal to the fraction ic(y) of the normal 
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density ratio e~ ay . The functions A , ix , p , a occurring in these equations are 
solutions of the first, second, fourth, and fifth of equations (12) with initial values 
( 1 , cot w , , ) at t = T , o) being the positive angle of fall. The function r is 
given by the equation (15). 

In constructing a range table it is sometimes desired to compute the first 
order correction to the range X due to the use of a slightly different value 
C + AC of the ballistic coefficient. The formula (18) will give this correction 
also, if the function n(y) is everywhere replaced by ic(y) — AC/C. For 
the equations (16) still hold provided that the expression (17) for F is multi- 
plied by the factor 

1 _ AC 

AC~ 1 C + '"' 
1 + C 

which gives, to first order terms, the value 



F=-U(F) e -^[l + K (2/)-^+...| 



C 



This is the same as before except that k ( y ) is replaced by k ( y ) — AC/C . 

The versatility of the formulas (13) and (14) of Section 2 is still further 
illustrated by the problem of determining the effect of the rotation of the 
earth upon a trajectory. Roever has set up the differential equations of 
motion in the form 

x" = — Fx' + 2o> ( y' cos sin a — z' sin 4> ) , 

(20) y" = — Fy' — g — 2u cos <j> ( z' cos a + x' sin a ) , 

z" = — Fz' + 2a> ( x' sin <f> + y' cos <f> cos a ) . 

In these formulas <f> is the latitude of the gun position, a is the azimuth of the 
plane of fire measured clockwise from the south, u = .0000729 is the angular 
velocity of the earth's rotation in radians per second, and F is the value (17) 
except that now w x (y) = w z (y) = n(y) = 0, since only the effect of the 
rotation of the earth apart from other disturbances is to be considered. The 
effects of all the disturbances could be evaluated simultaneously if desired 
but the result would be merely to add the terms about to be deduced to the 
expressions (18) and (19). In the equations (20) powers of a higher than the 
first have been neglected. 

The second member of the first of equations (20) is now the function 
f(t,x,y,z,x',y',z'), and the value of A/ in the formula (13) is the differ- 
ence / — / formed for the arguments t , x , y , z , x' , y' , z' of the initial trajec- 
tory (10), with similar values for Ag and Ah . It follows readily that 
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A/ = 2o>y' cos <t> sin a , 
Ag = — 2ux' cos <f> sin a , 
AA = 2w (a;' sin <{> + y' cos cos a) , 
and the formulas (13) and (14) give 

(21) 8X = 2w cos 4> sin a I (py' — ax')dt, 

Jo 

(22) 8z(T) = 2o> I t (x' sin <f> + y' cos $ cos a ) dt . 

Jo 

TAese are the first differentials of the range X and the deflection z(T) con- 
sidered as functions of the angular velocity co of the rotation of the earth. They 
give the first order effects of the rotation of the earth upon X and z(T) . 

The advantage, from the standpoint of the computer, of the formulas for 
the corrections derived in this section lies in the fact that only two integrations 
of systems of differential equations by methods of approximation have to be 
made in order to get all of the corrections. The first integration is of the 
system (9) in order to find the original trajectory (10), and the second is the 
integration of the first, second, fourth, and fifth of equations (12) for the 
functions (X, p., p, a) with initial values (1, cot w,0,0)atf=T. When 
this last set of functions has been determined the formulas (18), (19), (21), 
(22) give the values of the desired corrections by elementary operations of 
arithmetic, except in the case of the rotation of the earth for which two quad- 
ratures by methods of approximation are required. The details of the com- 
putation are described in the papers mentioned above* and need not be 
reproduced here. The methods in use when these results were obtained 
involved a separate approximate integration of a system of differential equa- 
tions for each type of correction, and this multiplicity of integrations is now 
replaced by two, the amount of computational work required being diminished 
by about three fourths. 

4. The differentials of the coordinates x, y , z for anti-aircraft 

CORRECTIONS 

The details of the computation of anti-aircraft corrections have been de- 
scribed by the writer in another paper, f Theoretically the steps used to 
obtain them are the following. 

A fundamental system of four solutions (Xj, p%, Pi, <?i) (i = 1 , • • • , 4) of 

the first, second, fourth, and fifth of equations (12) gives when substituted in 

equation (7) four linear relations 

♦Journal of the United States Artillery, vol. 51 (1919), p. 296. 

t Submitted for publication to the Journal of the United States Artillery. 
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(23) X i bx + m 5y + pi dx' + <?i dy' = m,i ( i = l , • • • , 4 ) . 

If the fundamental system chosen is the one whose elements reduce to the 
identity matrix at t = 0, the four values m t - from the second members of 
equation (7) are 

/»* /»* 

mi = I (pi A/ + o"i Ag)dt, m 2 = I (p 2 A/ + <r 2 Ag)dt, 
Jo Jo 

m z = Ax' + I (p 3 A/ + a 3 Ag)dt, 
Jo 

m i = Ay + I (P4A/ + <r 4 Ag)dt. 
Jo 

The differentials of x and y are then found by solving the four linear equations 
(23) for 5a; , by , bx' , by' , and these differentials are the first order corrections 
to x and y required at the time t on the trajectory. 

The effectiveness of this method of determining the differential corrections 
from the standpoint of the computer depends very much upon the plan of the 
computations and the arrangement of the work. Such a plan has already 
been described by the writer in the paper last referred to. It turns out that 
the fundamental system of solutions ( X, , in, p», <r») (i = 1 , • • • , 4) , and the 
cofactors of the determinant of the coefficients in the first members of equa- 
tions (23), can be determined by three approximate integrations plus ele- 
mentary arithmetical operations. The system of differential equations 
actually integrated is the adjoint of the first, second, fourth, and fifth equa- 
tions of the system (12) which has the form 

A' = C, 

B' = D, 



(24) 



C = ax' FB -F(l + x' 2 ^jC - x'y'F^D, 

D' = ay' FB - x' y' F^C - F (l + y' 2 ^JD . 

If (A, B,C,D) is & solution of this system, and ( X , ju , p , <r ) a solution of 
the original one, the relation 

^4X + B/x + Cp + D<r = const., 

is always true, and it follows readily that four solutions (Ai, Bi, C», Z)») 
(i = l,---,4)of the system (24), with initial values at t = forming the 
identity matrix, will have elements which are the cofactors of the corresponding 
elements of the determinant of the solutions (X,-, m, p», o-j) (i = 1 , • • • , 4) 
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of the original system. The matrix of solutions of the system (24) has the 
form 

10 
/2k\ ^2 B 2 C 2 D 2 

A3 B3 C3 D3 

Ai B4 Ci Di 

and it is clear from the form of the equations themselves that when the ele- 
ments B , C , D of one of the last three rows have been determined by inte- 
gration of the last three differential equations, the corresponding A can be 
determined by a quadrature. The elements of the solutions (X<, m, Pi, <r,) 
( i = 1 , • • • , 4 ) are the co-factors of the matrix (25) and they can be cal- 
culated by elementary operations of arithmetic. 

The differential of z is found by using in equation (7) two solutions ( v , t ) 
of the third and sixth of equations (12) with initial values ( vi , n ) = (1,0), 
(v2, t 2 ) = (0,1) at 2 = 0. These solutions are, respectively, 

x 

"1 = 1, ri = - -,, 

"2 = 0, 7 "2 = ;7> 

as may readily be established by elementary integrations similar to those 
used in deducing formula (15). Equation (7) then gives 



Sz — — , 



~- t sz' = - r -,Ahdt, 

Uz' = Az'o+ f ^Ahdt, 



Xn 
X 



which determine the differentials Sz and Sz' . 



